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The Van Vleck-Pauli-Morette fluctuation determinant is derived from the global group property 
of the time evolution amplitude in a continuous formulation of path integrals. 
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In the semiclassical limit, the quantum mechanical time evolution amplitude consists of an exponential of the 
classical action exp (i A) , multiplied by a fluctuation factor F containing the inverse square root of the functional 
fluctuation determinant of harmonic eigenmodes of the system. The standard derivations of F rely on the local 
group properties of the time evolution amplitude jjj. The initial historic paper of DeWitt-Morette (2) determined 
F by enforcing these properties for infinitesimal time slices of the amplitude, which are necessarily semiclassical by 
\ Dirac's observation The full fluctuation factor F for finite times was then composed by a limiting procedure from 
the F's of the time slices. The result was expressed as a square root of an ordinary matrix determinant, the Van 
<^"4 ' Vleck-Pauli-Morette determinant |2| . 

Cn| . Later, Gelfand and Yaglom Q related the fluctuation determinant to the solution of a second-order differential 
equation, again via time-slicing techniques. This solution can, of course, be related to the Van Vleck-Pauli-Morette 
determinant. 

In this note, we point out a more compact way of obtaining the fluctuation factor from the global, finite-time group 
property of the time evolution operator. Only the continuum formulation of path integrals is used. Our derivation 
. involves neither time-sliced actions nor differential equations. 
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O 1 II. SEMICLASSICAL APPROXIMATION 

o 

In Schrodinger theory, a point particle in a D-dimensional euclidean space has associated with it a Hamilton 
Oh, operator H(t), and a time evolution operator U(t), which determines the amplitude to go from a position x a at time 
-4— > ' t a to a position Xf, at time % by the matrix elements Q 
Ch 

(x b t b \x a t a ) = 0(t b - t a )(x b \U(t b , t a )\x a ). (1) 

a -1 

The Heaviside function 0(t b — t a ) ensures causality by the vanishing of the amplitude for times t b < t a . As elements 
of a one-parameter Lie group, the time evolution operators for different times satisfy the group multiplication law 

U(t b ,t a ) = U(t b ,t)U(t,t a ). (2) 

For matrix elements, this reads 

/OO 
d D x(x b \U{t b , t)\x)(x\U(t, t a )\x a ) (3) 
-OO 

so that the time evolution amplitudes satify the integral relation 

D 



(x b t b \x a t a ) =Y\ I dx 

■ 1 J — OO 



(X b t b \xt)(xt\x a t a )- (4) 



i=l 

The matrix elements in (fy) have a path integral representation 
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where 



(x b \U(t b ,t a )\x a ) = I Vx exp 



A{x) 



-A{x) 



dtL{x, x, t) 



(5) 



(6) 



is the action and L the Lagrangian of the system. The semiclassical approximation is defined rewriting the result of 
the path integration as a product 



(x b \U(t b ,t a )\x a ) w exp 



A(x£,, x ai i a ) 

a 



F^Xb) x ai ^5, £ a ) 



(7) 



where A(x b , x a ;t bl t a ) is the associated classical action, i.e., the action A[x] evaluated for the solution x c \(t) of the 
Euler-Lagrange classical equation of motion which extremizes A[x] with fixed endpoints at x a ,t a and x b ,t b . In the 
semiclassical approximation, the factor F contains no h and contains the fluctuation determinant arising from the 
quadratic fluctuations around the classical path. Its logarithm is the quantum-mechanical analog of the entropy of 
harmonic fluctuations in quantum statistical mechanics. Since the end points of the paths are fixed, x(t a ) — x a , 
x(t b ) — x b , the boundary conditions for the fluctuations Sx(t) are of the Dirichlet type: 5x(t a ) — 0, 5x(t b ) = 0. 
For a point particle moving in a time-dependent potential V(x,t), the Lagrangian reads 



L = ^-x 2 -V{x,t), 



and the fluctuation factor is 



F{x\)i x ai tfr, i a ) 



det D (d 2 /dt 2 ) 



-D 



(8) 



(9) 



y/ 2nih(t b -t a )/M D V det D [9 2 M 2 + VW(t)/M] 
where V (2) (f ) is a D x D derivative matrix collecting the second derivatives of the potential along the classical path 

d 2 



didj 



V(x,t) 



(10) 



X—X c l(t) 



where the indices i,j denote the vector components. The fluctuation determinants deto consist of the product of 
eigenvalues of the D x D differential operator for Dirichlet boundary conditions. 

The fluctuation determinant is most easily evaluated with the help of the Gelfand-Yaglom method H, and the 
result can be reexpressed in terms of the Van Vleck-Pauli-Morette determinant M 



F(x b ,x a :t b ,t a ) = /2 n ift\D/2 {detn[-d x i a d x ^A(x b , x a ;t b , t a )]j 



1/2 



(11) 



The minus sign inside the determinant makes the argument of the square root positive as long as the classical 
trajectories do not reach a turning point. The continuation to longer intervals can always be done with the help of 
Maslov indices [Q. 

There exists various ways of rewriting the right-hand side of Eq. ([ll]) . A convenient form is obtained by using the 
fact that the momentum at the final time is given by the derivative of the classical action: 



j uA[x b , x a ',t b ,t a ) 
PI = ~ 



dxl 



allowing to rewrite (|11|) as 



F(x bl x ai t b: i a ) 



(2TTlh) D / 2 



det D (-d K p 3 b ) 



1/2 



(12) 



(13) 



For the decomposition (0) of the matrix elements of the time evolution operator, the group property (||) takes the 
form 



,j-A(x b ,x a ;t b ,t a ) F ( , . \ 



D 

n 



e^^'^i^, x; t 6 , f) x e iA*( x , Xa ;t,t a ) F ( X! Xa . t> ta) (w) 



Here we have introduced superscripts Z/ and i? to emphasize the left and right positions of the action in the product. 
The configuration of the variables is illustrated in Fig. [I]. 




FIG. 1. The upper curve shows two classical paths running from x a to x and from x to x b . The intermediate point x has 
to be determined by the saddle point condition. The lower curve shows the direct classical path. 



The fluctuation factor ( ]11| ) has the property that the semiclassical approximation (|7j) satifies this equation providing 
that the intermediate x-integrals are evaluated in the saddle-point approximation. This operation is done explicitly 
as follows. We denote the extremum of the intermediate integration over x by x, and expand the integrand around x 
up to quadratic terms in 5x = x — x. Since the fluctuation factor contains no h, only the action in the exponent has 
to be expanded, and the semiclassical approximation to (pdh reads 



exp 



- A(x b: x a , tfr, t a ) 



— A L (xb, i; tb, t) + — A^(x, x a ; t, t a ) 



D 



F(x b ,x a ;t b ,t a ) = exp 

JJ d5x l exp 
-°° i=i 

The saddle point condition for x is 

|l(A^A-) = 0. 

Just as in Eq. (p"2|), the derivatives are equal to the momenta at the intermediate time t, 

■^-A R {x,x a ;t,t a )=p?(t), 

■J^-A L {x b ,x;t b ,t) = -pf{t), 
ax 1 

so that the saddle point condition ( |l6| ) implies the equality of the intermediate momenta: 

(*)=*?(*)• 

Our proof will be straightforward for a general Lagrangian which is at most quadratic in the velocitie: 
L(x,x;t) = ^x l g lj (x,t)x J +x l a l (x,t) ~ V(x,t), (g^ = g ri .) 



(15) 

(16) 

(17) 
(18) 

(19) 

(20) 



The kinetic metric gij(x,t) is also known as the Hessian, whose determinant is assumed to be nonzero to have a 
nondegenerate quantum system. 
The canonical momenta are 



dL _ . 

Pi — -qT — 9ij Xj + a,; 



and the Hamiltonian is the Legendre transform of L(x, x; t): 

H(t) =p t x l -L = -\pi(x,t) - ai{x,t)]g ij (x,t)\pj(x,t) - aj(x,t)] + V{x,t), 



(21) 



(22) 



where g 1J {x,t) is the inverse matrix of the Hessian gij(x,t). 
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The Euler-Lagrange equations of motion following from (|20|) are 

d dL dL 



dt d&i dx, 



= 0. 



(23) 



They are second-order differential equations in time. Due to this property, the condition ( |l9| ) at the junction between 
the left and right paths ensures that the saddle point x is located on the the classical trajectory x c \{t) running all the 
way from x a to xt,. The upper curve of Fig. [I] coincide then with the lower one. Moreover, the sum of the actions is 
equal to the total classical action along this combined path: 



A^{%bi 5>] tb, t) + A (x, X a ] t, t a ) — A(xb, X a ] tb, t a ). 



(24) 

Performing the Gaussian integral in Eq. (|l^), we therefore find the semiclassical consequence of the group property 
(pTh for of the fluctuation factor: 



F(x b ,x a ;tb,t a ) = F(x b ,x;tb,t)F(x,x a ;t,t a )(2TTih) D ^ 2 < det 



d 2 (A L +A R ) 



dx l dx : > 



-1/2 



(25) 



where we have omitted the tilde on top of the intermediate position x on the classical path. This equation is an 
algebraic version of the eikonal equation in Schrodinger theory. 

It is straightforward to verify that the Van Vleck-Pauli-Morette fluctuation factor ( pd| ) satisfies (|2^) . We shall prove 
this using the equivalent form (|li|). Inserting this into (|5j), and using ([Tt| ) and (|l^), we have to show that 



det; 



dpb 
dx a 



det D (- dpZ/dxl^ det D (- dp R /dx a \ x ) 
det D (dp R /dx\ Xa - dp L /dx\ x J 



(26) 



where we have ignored vector indices, for simplicity, and emphasized the variable kept constant in the partial differ- 
entiations. We also used p L — pb- The proof follows from the chain rule for the Jacobians. Taking the left-hand side 
of Eq. (|2^) to the right-hand side, we must verify that 



f 



det 



D 



dp R /dx a \J detD (dx a /dx\. 



det D (dp R /dx\ Xa - dp L /dx\ Xb ) 
The saddle-point condition (|l9|) is, in a more explicit notation, 

p L [x b ,x(x a ,x b )} = p R [x a ,x(x a ,x b )}. 

This equality allows to derive 



d P L 




dp R 


_ dp R [x a ,x(x a ,Xb)} 




dp R 


dx a 




dp R 


dx 


XI, 


dx 


9X 




8x a 


x dx 


H 

Xb 


dx 



(27) 



(28) 



(29) 



Inserting this equation into the denominator of (ETj) proves that (j26|) is indeed satisfied. 



III. GLOBAL DERIVATION OF FLUCTUATION FACTOR 



We are now prepared for the essential part of this paper, in which we derive the Van Vleck-Pauli-Morette formula 
( |i"l| ) from the semiclassical group property ( p5| ) of the fluctuation factor. We proceed in two steps: first we move the 
intermediate time t infinitesimally close to the initial time t a . This time is called t a + . The corresponding intermediate 
position x will then lie at a point illustrated in Fig. ||. For this configuration, the fluctuation factor 

( pB| ) reads more explicitly 




FIG. 2. Classical path from x a to an intermediate position x — x a + very close to x a , followed by a classical path from x a + 

tO Xb- 
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F(xb, x a ;t b , t a ) = F(x b , x a + ; i b , t a +)F(x a + , x a ; t a + , t a ){2mh) 



D/2 



x < det 



D 



d 2 A L (x b ,x a +- 1 t b ,t a +) d 2 A R (x a +,x a ;t a +,t a ) 



dx l dxi 



dx l dxi 



-1/2 



(30) 



In the limit t a + — > t a we now extract the behavior of F(x a + , x a ', t a + ,t a ). Intuitively, this should be determined by the 
kinetic term of the action only, since the potential has no time to become active. Let us see how this comes about. 
First we use the equations for the momenta (|l7|) and (|l8|), express these in terms of the derivatives of the Lagranian 
with respect to the velocities via (|2l|), and derive the relations 



o 2 l 



d 2 L d 



d 2 



dx l dx3 



A L (x bl x]t b ,t) 



_d 

dx^ 



Pi(t) 



dxidij dxkdij dx 
d 2 L d 2 L d 



dxidij dikdij dx 



■x L k (t). 



(31) 
(32) 



On the right-hand side we have taken into account that the arguments x and x of the Lagrangian are classical 
trajectories fixed by their end points, i.e., x(t) = x c \(x,x a ;t) in Eq. ( |3l| ) and x(t) — x c \(x b ,x\t) in Eq. (|32"|). The 
derivatives with respect to the end points produce therefore an extra term coming from the velocity dependence of 
L(x,x;t). Then we use the Lagrangian (EOT) once more to express 



dL 1 

-q^T = -± l [d k g l j(x 1 t)}x j + x i d k a l (x,t) - d k V(x,t) 



d 2 L 

diidxj 

d 2 L 
dxidxj 



[djgik(x,t)]x k +dja,i(x,t) 



9ij{x,t) 



such that the brackets in Eq. (|30| ) lead, via (|3l| ) and (32), to 
d 2 A R {x,x a ;t,t a ) d 2 A L {x b ,x;t b ,t) 



dx l dxi 



dx l dxi 



(33) 
(34) 
(35) 

(36) 



d d 

= {[dj9ik(x,t)]x k (t) +d j a i (x,t)} R - {[djg ik (x,t)]x k (t) +d J a l (x,t)} L + g lk (x,t) R —x R (t) - g ik (x,t) L —x k (t). 

Then we use the fact that x(t) and x(t) are continuous at the junction between the left and right paths, such that we 
can collect the terms on the right-hand side to 



d 2 A R {x,x a ;t,t a ) d 2 A L (x bl x]t b ,t) 



dx l dxi 



dx i dxi 



9ik(x,t) 



i^-x R (t) - -^-x k (t) 

dxi k w dx^ 



(37) 



If we now take the limit t 



tn 



the contribution from the path R to the derivatives inside the brackets becomes 



much larger than that of the path L. Indeed, the associated short classical path is 

Xk(t a +) - Xkita) 



tn 



implying a very large derivative 



ta 



t„ 



ta 



(38) 



(39) 



Inserting this dominant contribution into Eq. (p7[), this further into (|3C|), and factorizing out the approximately equal 
unknown fluctuation factors F(x b , x a ;t b , t a ) « F(x b , x a + ; t b , t a +), we obtain the fluctuation factor for the infinitesimal 
time interval: 



F(x a + , x a , t a + , i a ) 



1 



(2inh) D / 2 



det 



D 



gij \Xa , t a ) 



tn 



t a 



1/2 



(40) 



This is the well-known free-particle result, as anticipated. It will be used twice: first to obtain Eq. (|41j), and later to 
fix the sign of the solution in Eq. (p3h. 
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We turn now to the second step in the derivation. It is based on the observation that the group property ( pq ) is 
not only valid for causal time configurations tb > t > t a , but also for acausal time configurations t > tb > t a . The 
caus ality is only a property of the time evolution amplitude ([!]) , not of the matrix elements (]?]) , such that also jl^ ) 
and ( |25| ) are valid for tj, > t > t a and t > tb > t a . This means that we can also bring the time tb close to t a , leaving t 
much larger than these two adjacent times, as indicated in Fig. H. 




FIG. 3. Classical path from x a to x, followed by a classical path backward in time from x to Xb along the same path. 



In this limit tb — > t a + , Eq. (£5|) reads 

F(x a +,x a ;t a +,t a ) = F(x a +,x;t a +,t)F(x,x a ;t,t a )(2Trih) D/ ' 2 \ dct D 



d 2 (A L +A R ) 



dx % dxi 



-1/2 



where A L and A R are now abbreviations for A L = A(x a + , x; t a + , t), A R = A(x,x a ;t,t a ). Using (B 
rewritten as 



(41) 

this can be 



i 



(2mh) D 



det / 



9ij [p^a i ta ) 



t a 



1/2 



det; 



d 2 (A R + A L ) 
dx l dxi 



1/2 



Let us study the behavior of the brackets in the last determinant for t a + close to t a . It reads more explicitly 

d 2 A R (x,x a ;t,t a ) d 2 A L (x a+ ,x;t a+ ,t) 



dx i dx> 



dx l dxi 



(42) 



(43) 



where x a + is very close to x a . As the limit t a + — * t is reached, the two paths coincide, and have the same classical 
action, except for a negative relative sign, since the corresponding paths have the opposite direction in time. Thus 
for t a + — t a , the sum in ( [43| ) vanishes. For small t a + — t, we perform a Taylor expansion of the second term around 
the first and have, omitting the now superfluous distinction between L and R, and using double primes to abbreviate 
the second derivatives d /dxidxj, 

A" ix a + ? x i t a + j t) ~ A" (x a , X\ t a , t) + — A (x a , X] t a , t) (t a + t a ) + 



d 



( X a+ — X a) 



(44) 



Inserting here x^+ — x k ~ (t a + — t a ) xk {t a ), we may replace ( (43| ) by 

d d 2 A(x a ,x;t a ,t)] \ d d 2 A{x a ,x]t a ,t) 



dt a dx % dx3 



dx% dx l dxi 



i k {ta) > (t a +-t a ), 



(45) 



and 



becomes 



F(x,x a ;t,t a )F(x a ,x;t a7 t) = 1 {detg Iffy i x a, t a )}} 1/2 



x det 



D 



d 2 dA(x a ,x;t a ,t) 



dx x dx 3 



dt a 



d 2 dA(x a ,x;t a ,t) 



dx l dxi 



dx„ 



x k (t 



■>}) 



1/2 



(46) 



Note that the derivative d 2 /dx l dx^ does not act on x k (t a ), so that the total argument in the last determinant of ( |4^ ) 
is not the double prime of the total derivative of the action with respect to time t a , in which case it could have been 
simplified to (d /dx l dx 1 )dA{x a , x; t a , t)/dt a = (d 2 /dx l dx-')L(x a , x(t a ), t a ). Since this is not the case, we can only do 
a partial simplification using the Hamilton- Jacobi equation 



d 

— A(x a , x; t a , t) + E(x a ,x; t a , t) = 0, 

dt a 



(47) 



7 



where E(x a , x; t a , t) is the energy at the time t a for the classical trajectory x(t a ) = x c \{x a , x; t a , t) running backwards 
from x at t to x a at t a . It is the value of the Hamiltonian H(t a ) of Eq. ( |22|) evaluated for this trajectory. If 
p(ia) = Pc\{x a ,x;t a ,t) denotes the associated trajectory in momentum space, the energy is given by 



E(x a , x; t a , t) — H(t)\ x ( t j =Xcl ( Xa X . ta t -) p ( t } =Pcl ( Xa X . ta ^, 

= ^[pi(x a ,t a ) - a,i(x a ,t a )]g v (x a ,ta)\pj(x a ,ta) - a,j(x a ,t a )] + V(x a ,t a ). 
With this, Eq. @ allows to rewrite @ as 

{det D [gij{x a ,t a )]} 1/2 



(48) 



F{jc^ x a , tdjF "(x , x, £ a , — 



^det D | 



{2Trih) D 
d 2 



E {x a ■ x . t (i ^ £) 



<9 2 dA(x a ,x;t a ,t) 



dx l dxi 



dx k „ 



x k {t 



1/2 



(49) 



At this point we observe that for a purely harmonic Lagrangian, which is at most quadratic in the velocities and 
positions, the functions gij(x,t) in the general expression ( p0| ) are position-independent, the vector potential cii(x,t) 
is at most linear in x, and the scalar potential V(x,t) has the general form V = Xiflij(t)xj /2. Then the classical 
action is at most quadratic in the end points. This implies a vanishing second term in the brackets of the second 
determinant in (E9h. Then, we have 



F{x, x ai t, t a ^f 



{det D Mi a )]} 1/2 



det/ 



dx l dxi 



E i^Xa i X>\ ta t) 



1/2 



for quadratic actions (^lltiK)^ 

Note the sign change of the second derivative of E{x a -> x; t ai t). This is caused by the replacement 

F(xj x ai t a ^F{x aj t aj £) > i F^x, x a , t 7 t a ) . 



(50) 



(51) 



The reason for the factor i D lies in the Fresnel nature of the path integral over the fluctuations. The exponent is 
the second functional derivative of the action with a factor i. Assuming stable orbits, the factor of i is positive or 
negative, depending on the time direction of the path. This sign change implies that the Fresnel integrals are related 
by 



F(x,x a ;t,t a )Vi — F(x a ,x;t a ,t)/Vi ■ 



We can now take the square root of ( pOf ) and obtain 

{det D [ 9lJ (t a )}} 1/4 



F{xhi %a'i tb-i to) 



for quadratic actions (27ri/i)^/^ 



det i 



dx l b dx 3 b 



E {x a , Xb , t a , ib) 



1/4 



(52) 



(53) 



The sign of this square root is fixed by the fact that in the limit of short intervals tb — t a , the fluctuation factor has 
to reduce to the free-particle result (ffol). 

Note that in the semiclassical limit the fluctuations are always harmonic. For a vanishing vector potential a,i(x,t) 
in Eq. (^o|), these would be driven by a time-dependent frequency matrix 



= g ik {x c i(t),t)d k djV{x cl (t),t). 



(54) 



This harmonic property does not, however, allow us to use formula ( |53|) for the fluctuation factor, since the frequency 
matrix depends on the end points via the classical solution x c \(t) of the equations of motion, so that the full formula 
( [49| ) must be used, which we now investigate in detail. 
We must evaluate 



d 2 



dx l dx : > 



E(x a ,x;t a ,t) + 



d 2 dk(x a ,x;t a ,t) 



dx l dxJ 



dx^ 



with dA(x a ,x;t a ,t)/dx a — Pk(t a )- Since dx a /dx 3 — 0, expression (B5h can be rewritten as 



1 



d 2 



Xk(t a )xi(t a ) + x k (t a )g k i(x a ,t a ) 



d 2 



dx l dxi 



xi(t a ). 



(55) 



(56) 



Using the symmetry gki{x a , t a ) = 9ik(x a ,t a ), this becomes 



, d , . . d . , 

9kl{Xa,ta)-Q — X k {t a )—Xl(t a ). 



Now the determinant of a product of matrices factorizes into a product of determinants, and Eq. ([49| ) becomes 



F{x, x a: t 7 i a ) 



:2 _ detp [gij{x a ,tg)] 
(2irih) D 



det 



D 



d . 

dx iXj[ta > 



where the phase factor ( |5l| ) has been taken into account. Using the relation 

d 2 d d 

■^—^—A(x,X a ;t,t a ) = --Q-jP^a) = -gik(x a ,t a )—X k (t a ) 



we can finally rewrite (pq) as 



F{x, x a , t, £ a ) 



from which it is straightforward to obtain 



■ det; 



d 2 



dxl dxi 



A(#, x a 5 £ a ) 



(27T^/ 2 



det 



D 



"1 ^ 1/2 



dx l a dx J b 



A(x^, x a , t a ) 



(57) 



(58) 



(59) 



(60) 



(61) 



This is the Van Vleck-Pauli-Morette formula @. 

Our derivation has ignored zero modes in the intermediate integration, which may be treated in the standard way 



IV. APPLICATIONS 



Here we shall apply our formula to three systems of point particles: 
A. the free point particle with a mass matrix, with a Lagrangian 

L —XindijXj) 



(62) 



B. the harmonic oscillator space with a time-dependent frequency matrix u>ij (t) , a mass matrix and a Lagrangian 



L = - {xiMijXj - Xi[Mu) 2 (t)]ijXj} 



(63) 



C. an ordinary particle in a constant magnetic field perpendicular to the plan spanned by x\ and x%, with a 
Lagrangian 



1 ° 
L = L= -M2J 



X; B±2X\. 



(64) 



In each cases, the boundary conditions are x(t a ) — x a ,x{%) = £&• 



A. Free Particle 



The free particle case is particularily simple. The Lagrangian (p3) implies the equations of motion 



M^x 3 =0. 



(65) 
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Since the matrix M is symmetic, it can be diagonalized by a similarity transformation with an orthogonal matrix S= 
S _T . Let M d = S _1 MS be the resulting diagonal mass matrix. The normal modes of the motion are y(t) = S -1 x(t). 
The latter satisfy 



dt 



and have the time dependence 



1 



Ui(t) = —— [y l a(t b -t)+yt(t- t a )] . 



h — ta 



The associated classical Hamiltonian is 



and the trajectories have the energy 



H = -x i M lj x 2 = -ViMfjitj, 



E(x a ,x b ;t a ,t b ) = - 



k=l 



{h - t a f 



Using the relation 



we deduce 



d _ dy^ d _ ! d 



dx\ dx\ dy b 



= Q~ 



O 2 



dx\dx > h 



E{x ai x b \ t a: t b j 



Mi. 



{t b -ta) 2 ' 

Inserting this into Eq. (53), we obtain the well-known fluctuation factor 

(det D M) 1/2 



[2mh(t b - t a )} 



D/2 ' 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



B. Harmonic Oscillator with Time-Dependent Frequency 

The case of the harmonic oscillator with a time-dependent frequency is slightly more involved. One cannot solve 
the equations of motion to get the solution in a closed form. We will however give a formal solution, showing how 
the well-known result can be recovered when the frequency is time independent. The equations of motion associated 
with ([33]) are 

d 2 

^Si + <4(t)ay=0. (73) 

With the help of two matrices A and B, the solution can be decomposed as x = A-x a + B.tj,. Since x a and x b are 
independent, each of the matrices satisfies a same equation as (f73"[): 

^A + a; 2 (i)A = 0, (74) 
^B + U 2 (t)B,0. (75) 



The boundary conditions are 



A ij (t a ) = 8 ij ,A ij (t b )=0, (76) 
B ij (t a ) = 0,B ij (t b )=6 ij . (77) 
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Formula (53) contains a double derivative with respect to the end point Xb- For this reason, it will depend only on 
the part of the solution with the matrix B. This simplifies the evaluation of the Hamiltonian and, taking derivatives 
with respect to the end point Xb, we have 



d 2 



dx l b dx 3 b 



E{x ai Xbi t ai tb) 



dt 



B(i) 



M 



B(i) 



B(t )MB(f ) 



(78) 



where the last equality defines B as the time derivative of the matrix B. Using this relation in Eq. (|53|), the fluctuation 
factor is given by 



(det D M 



1/2 



[2mh(t b -t a )] D/2 



det D B(t a ) 



1/2 



(79) 



which requires to solve Eq. ([rH) with the associated boundary conditions ([771). A formal solution can be obtained in 
the following way. Integrating twice (fr5|) , using ^7%, leads to 



B(i) 



dsB(t a ) 



ds / uj 2 {y)B{y)dy 



which can be iterated to lead to a Neumann series 



B(t) 



ds < 1 



0J 2 {y)dy ds'+ LJ 2 {y)dy ds' u 2 {y')dy>\ ds" 



B(f ). 



(80) 



(81) 



Using a first order differential formalism, this expansion can be given a compact notation. This comes from the 
fact that the solution of (fTa) can be written as 



B(t) 
B(t) 



T <^ cosh 



1 

-LJ 2 (s) 



f/.S 



B(t„) 
B(t 8 ) 



T < sinh 



1 

-w 2 (s) 



rf.s 



B(t„) 
B(i„) 



(82) 



where the hyperbolic functions are defined through their Taylor expansion and where the symbol T implies a time 
ordering operation. Using the boundary conditions, we end up with 



B(i) 




= T < sinh 



to 



1 

-u?(s) 





B(t a ) 



(83) 



where we have multiplied from the right by an appropriate matrix in order to single out the upper left component. 
We can extract B(t a ) from this relation using the boundary condition B(tb) = 1. This gives 



B(f„) 



ds 



(84) 



which is indeed equivalent to iterative solution (^l|). Using this in (|7^) provides then us with a formal solution for 
the fluctuation factor. 

The case of a time independent frequency matrix is obtained directly from the formal series ( |84| ) : the time ordering 
operator disappears, the time integration is trivial and the series easily evaluated: 



T I sinh 



1 

-LJ 2 (s) 



ds 



sin [u)(t b - t a )] 



u>~ 
-v 



As stipulated in (B4J), we need only the upper- right component. We then end up with 



F(Xb, Xa> tb, ta) 



(detzjM) 



1/2 



det/ja; 



1/2 



a;=co„st. [2mh(t b - t a )] D/2 I dct D sm[u(t b - Q] 
Denoting by u 2 (no summation over i) the normal modes, this equation can also be written as 

,1/2 D r n 1/2 



(det^M) 1 



OJ— const. 



n 



[2mh(t b - t a )] D/2 fJi I sin[u>i(t b ~ t a )} 



(85) 



(86) 



(87) 
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C. Particle in Constant Magnetic Field 



Here the calculation looks somewhat more complicated, although it is still trivial. The equations of motion associ- 
ated with the Lagrangian (fil) are 



d 2 d 
d 2 d 

_ X2 + aJ _ Xl = 0> 



(88) 
(89) 
(90) 



The index j will be limited to j — 3, ■ • ■ , D throughout in this section. The frequency lj is the Larmor frequency 
u> = eB/(cM), where e is the electron charge and c the speed of light. The classical trajectories are 



Xi 



■r-2 



sin u(tb - t a ) 
1 

sinw(t b - t a ) 
1 



[{x\ - x\) sinw(t - t a ) + [x\ - Xq) smuj(t b - t)] + xj, 
[(xl - xl) sincj(t - t a ) + (x 2 a - Xq) smu>(t b - t)] + x%, 



H — ta 



Xi(t b - t) + xl(t - t a ) 



where and Xq are determined from pH ) and ( p9[ ) as 

' Li+4) + (^-^)oot^^ 



1 



(^+^)-(^-^)cot^^ 



For the classical Hamiltonian we obtain the only non-vanishing contributions (no summation over j) 

MU! 2 



d 2 d 2 
■—^-—^E(x a ,Xb;t a ,tb) — q x 2q x 2 E(x a ,x ] t a , t ) 



4sin 2 [a;(i fc -t a )/2]' 



i) 2 



dx J b dx J b 



EyX a: £ a , 



M 



from which it is trivial to find 

d 2 



dct 



D 



dx\dx^ 



- E {Xa , Xb, t ai tb) — 



M 


D-2 | 


[(tb-t a ) 2 _ 





Moj 2 



.4 sin 2 [w(t 6 - 1„)/2] J ' 



such that Eq. @ yields §J 



-D 



E(Xbi *^a; ^6? ^a) 



M 



u{t b -t a )/2 

2mh(tb - t a ) sin [u>(t b - t a )/2] ' 



(91) 
(92) 
(93) 

(94) 
(95) 

(96) 
(97) 

(98) 



(99) 



D. Particle in Arbitrary One-Dimensional Potential 



For a particle moving in an arbitrary time-dependent potential in one dimension, it is possible to construct an 
explicit solution to the general relation (^5|). The matrix in the determinant on the right-hand side is 



(100) 



Using the equation of motion 
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dt 2 



x + d x V{x,t) = 



(101) 



and taking the derivative with respect to time we obtain 



— x + V (x,t)x = 0, 
and a similar equation for a derivative with respect to any other parameter A: 



(102) 



— dxx + V {x,t)d x x = 0. 



(103) 



Hence d\x can be expressed as a linear combination of two fundamantal solutions of ( |102| ). One of them is the time 
derivative of the classical trajectory, d\x^ = x. The other can be obtain from the D'Alembert construction PJ and 

is d\x^ = x dt/x 2 . Combining these into solutions satisfying the boundary conditions d x x{t a ) = 0, d x x(t D ) = 
and d x x(t) = 1, we then obtain 



dxX 1 



drX 1 



x R jjdt/x 2 
i(t) fldi/i 2 ' 

± L J* b dt/x 2 
i(t) J*» dt/x 2 ' 



from which we deduce 



d x x R {t)-d x x L {t) 



ft dt/x* 



(j t [dt/x 2 ) (j t tb dt/a 
Inserting this result in (p5[), we have in the limit t — > t a + 

-1/2 



±\t a ) 



dt 

72 



= V 2nihF{x a + ,x a ;t a +,t a ). 



Taking the limit tb —> t a at an arbitrary t in (|25|), and inserting (107) , we obtain 

-.-1/2 



x 2 (t a )x 2 (t) 
and thus a fluctuation factor 



dt 



dt 



= (2ltiK)F(x a ,X] t a ,t)F(x, X a ;t, t a ). 



F(Xb, t 0l t a ) 



x(t a )x{t b ) 



dt 



-1/2 



(104) 
(105) 

(106) 



(107) 



(108) 



(109) 



which agrees, of course, with formula (O 



V. CONCLUSION 

We have shown that the Van Vleck-Pauli-Morette determinant in the fluctuation factor ( |TT| ) can be obtained 
directly from the group property of the time evolution operator (Q) and the semiclassical expansion (^). In addition, 
we have derived a formula which allows us to find the fluctuation factor form the classical Hamiltonian function if the 
Lagrangian is at harmonic velocities and positions. 
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